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SUMMARY 




An experimental investigation is carried out to determine the 
effect of an axially symmetric initial imperfection on the buckling 
load of a circular cylindrical shell under axial compression loading. 
The imperfection studied has the form of a half sine wave in the length 
direction. The shells are formed by a copper electroforming process. 
The theoretical buckling loads are also computed for the imperfect 
shells, and show the same trend as the experimental study. 


INTRODUCTION 


The large discrepancy between the classical buckling load of a 
circular cylindrical shell under axial compression and the experimental 
values of the buckling load is a well known but little understood phenom- 
enon (Ref. 1). Several reasons for this discrepancy have been proposed 
but have yet to be substantiated by adequate experimental investigations. 
Donnell and Wan (Refs. 2 and 3) have shown that imperfections in 
geometrical shape have an important effect but the exact degree of this 
effect for various types of imperfections is still to be determined. The 
idea of a jump of the shell from the prebuckled state at a load lower than 
the classical load has been proposed by von Kkrmkn and Tsien (Refs. 4 
and 5) but the exact mechanism of the jump and conditions under which it 
is possible are as yet unexplored. 

The purpose of this investigation was to study the effect of a 
specific type of initial imperfection under very carefully controlled 
conditions. In order to accomplish this the following steps have been 
carried out: 

1. development of a method of making cylindrical shells 
which were as nearly perfect as possible: 



2. establishment of testing techniques in which the stress 
distribution was controllable and measurable: 

3. determination of the buckling stress of these nearly 
perfect shells: 

4. addition of the initial imperfections to be studied and 
the determination of the buckling stress. 

The investigation was conducted at the Graduate Aeronautical 
Laboratories of the California Institute of Technology under the 
sponsorship and financial assistance of the National Aeronautics and 
Space Administration, Research Grant NsG-18-59. 
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EXPERIMENTAL INVESTIGATION 
Fabrication of the Test Specimens 


The method of manufacturing of the cylindrical shells is not new 
with GALCIT (see reference 6) but a number of refinements were 
required in order to obtain the uniformity and accuracy desired. Basi- 
cally it consisted of plating a copper shell on an accurately machined 
wax mandrel and then melting the mandrel out of the shell. For a shell 
8 inches in diameter and 10 inches long, a steel cylinder 7 inches in 
diameter and 13 inches long was used as a core (this provided a means 
of water cooling the wax to harden it). On this core was cast a layer of 
wax consisting of a two to one mixture of refined paraffin and Mobile 
Cerese Wax 2305. This was cooled, machined to the proper shape, and 
sprayed with silver paint thinned with toluene. Figure 1 shows a mandrel 
and finished wax form. The plating was done in a Cupric Fluoborate, 
Cu(BF 4 ) 2 > bath with a 15 inch diameter copper anode which was bagged 
in fine mesh Dynel fabric. During plating the mandrel was rotated and 
the bath was given additional agitation by forcing air through it. Using 
voltages less than 10 volts and current densities up to 55 amperes per 
square foot, the plating time was approximately 20 minutes per 0. 001 
inches of plate. 

After plating, the ends of the cylinder were removed (some 
thickness increase occurred just at the ends), the wax was carefully 
melted out, and the residual wax and silver paint was removed with 
benzene. The shell thickness used in this series of tests was approxi- 
mately 0. 0045 inches. The thickness variation in any shell could be 
held to less than 3 per cent. In the radial direction, the desired radius 
could be held to + t/2 or approximately + 0. 002 5 inches. 
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Tests to determine the characteristics of the plated copper were 
carried out in uniaxial tension. This was done by utilizing long strips 
of the copper that were soldered into one-eighth inch thick plates that 
were in turn clamped into the jaws of an Instron testing machine. The 
strips had length to width ratios greater than 15. The head displacement 
of the testing machine was used as the measure of strain and the load 
was read from the Instron load cell. A typical stres s- strain curve is 
shown in Figure 2. The value of Young's modulus used in this work was 

13. 0 x 10 6 lb/in 2 . 

A determination of Poisson's ratio was not attempted since its 
influence in the reduction of the buckling data is of secondary importance. 
A value of 0. 30 was used for this purpose. 


Test Procedure 


The buckling tests were carried out in the controlled displacement 
testing machine shown in Figure 3. This machine was designed to be 
rigid in comparison with the test specimen and capable of subjecting the 
test cylinder to very small increments of end displacement. The relative 
displacement of the two end plates of the testing machine was controlled 
by three screws. One complete turn of the screws gave a displacement 
of 0. 025 inches. The screws could be operated independently to give the 
proper load distribution on the shell and then simultaneously to increase 
the load to the buckling point. The springs shown in the figure were used 
to preload the testing machine when mounting the test specimen in the 
machine and securing it to the end plate of the testing machine. The 
testing was carried out when the machine was in the position shown in 
Figure 3. The end plate with the gear drive rested on pins and the 
opposite end rested on a set of rollers. 

The load distribution was monitored and total load was obtained 
on the load measuring cylinder shown in Figures 4 and 5. This consisted 
of a seamless brass cylinder which was 0. 0107 inches thick, 2. 5 inches 
long and 8. 000 inches in diameter. Twenty-four foil-type strain gauges 
were mounted on the cylinder equally spaced around the circumference. 
The ones on the inside were directly opposite to those on the outside. 

In order to see if the strain being measured in the load measuring 
cylinder corresponded to the strain in the test shell, the following test 
was carried out. A test cylinder was instrumented with strain gauges 
that were located at every other circumferential position as those on the 
load measuring cylinder. The test shell was mounted in the load 
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measuring cylinder and both sets of gauges read as the shell was 
eccentrically loaded. Two typical readings showing the correspondence 
of the measured strains are shown in Figure 6, 

The test cylinder was cast with a low melting temperature alloy 
into an end ring and the other end was secured to the load measuring 
cylinder in the same manner. After this operation was completed the 
shell was measured to determine the initial imperfection. 

In this series of tests the initial imperfection that was measured 
was the deviation of the generators of the shell from a straight line. 

This measurement was carried out at nine stations around the circum- 
ference. The measurements were made with a reluctance type pick up 
which has a sensitivity of about 2 5 volts per inch. The noise level and 

_4 

drift were such that deflection of 10 inches could be accurately read 
without making contact with the measured surface 

The measurements were carried out by mounting the sensing 
head of the pick up in a slide which traveled on a guide. The guide was 
attached to the end rings of the shell to be measured and readings were 
taken at 32 stations along the generator. The pick up and guide are 
shown in Figure 7. 

After the initial imperfections were measured, the end ring was 
then secured to the end plate of the testing machine with a thin layer of 
Devcon. After hardening of the Devcon was complete, the cylinder was 
ready for testing. 

The buckling test was carried out in the following manner. The 
cylinder was initially loaded to about one -third of the expected buckling 
load, and the circumferential load distribution was made as uniform as 
possible by adjusting the three screws of the testing machine. The load 
was then gradually increased in small increments by turning the three 
screws simultaneously. After each increase the load distribution was 
adjusted again. This was carried out up to about two-thirds of the 
expected buckling load. After this point the load distribution was not 
adjusted so as to prevent buckling occurring during one of the adjust- 
ments. The load was increased in small increments and the strain gauges 
monitored until buckling occurred. 
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Test Results 


Two types of axially symmetric imperfections were tested in 
this series of tests. Thirty- seven shells with an initial deformation 
of the form AR = a sin irx/ L were buckled and three with an imper- 
fection of constant curvature given by AR = a Q £2(2x/L) - (2x/L) z J 
were also tested. Table I shows the intended imperfection of the shells 
and Figures 8 through 16 give the measured imperfection on several of 
the shells. These figures show the shell generator at nine locations 
equally spaced around the circumference. Table I also shows the wall 
thickness of the shell tested. All shells had a base diameter of 8 inches 
and a length of 10 inches. 

Table I shows that some of the shells had an initial buckling. 

This initial buckling consisted of the formation of one wave on the 
surface of the shell. This caused the load distribution to fall off in the 
neighborhood of the one wave without appreciably affecting the distribu- 
tion over the rest of the circumference of the shell. After the initial 
buckling, the load was increased until general collapse occurred without 
attempting to alter the load distribution. General collapse occurred in 
the same manner as for the shells that did not have an initial buckling. 

General collapse consisted of a snap-through which is character- 
istic of this type of testing. In all but a few cases the post-buckled state 
consisted of 2 to 3 rows of buckles that extended completely around the 
circumference. The number of circumferential waves is noted in Table 
I. All of the shells that had an initial imperfection amplitude greater 
than 0. 010 inches buckled at one end or the other. This buckled state 
consisted of 2 to 3 rows of buckles that started quite close to one end 
and extended about one-third of the way up the shell. There did not 
seem to be a preferred end for this buckling to take place. The other 
shells with small positive imperfection, negative imperfection, and no 
intended imperfection buckled over the middle third of the shell. Again 
this consisted of 2 to 3 rows that extended completely around the 
circumference. 

Table I also gives the maximum variation in load distribution 
near buckling. The average maximum variation in load distribution 
for all the shells tested was 18. 6 per cent. In most cases the average 
stress was nearer the maximum stress rather than the minimum. 
Figures 17 and 18 show the load distribution on two shells as the load 
was increased. The initial buckling is easily detected in Figure 18. 
Figures 19 and 20 give the load distribution on several of the other 
shells at the last reading before buckling. 
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The values of K for this series of tests are given in Table I 

and the values of cr / cr are plotted versus a / 1 in Figure 21. 

cr ci o 

These are compared with the Kanemitsu and Nojima value which for 
the L/R and R/t corresponding to these specimens gives K = 0. 1 7. 
The theoretical curve shown in the figure is the result of the analysis 
given in the Appendix. The following important features can be seen: 

1. That with proper care in manufacturing and testing, 
values of the buckling stresses can be obtained 
which are much higher than those usually found. 

2. That, for the displacement forms tested, small 
departures from initial straightness lower the 
buckling stress and that the effect for inward 
displacements is greater than that for outward 
displacements. 

3. That if the outward displacements are increased 
the value of the buckling stress again rises until it 
reaches essentially the same value as that for the 
initially straight cylinder. 

4. That the constant curvature and sine wave shape 
give essentially the same values of cr for the 
larger values of a Q /t. 

The drop in the value of cr for 0 < a Q /t <10 and the return 
to the "perfect" value for a /t>2(fcannot be explained at this time. It 
is felt that some of the scatter in the experimental results can be 
accounted for by the variation of stress around the circumference of the 
shells. Also the lowering of the value of u below the classical value 
may be in part due to the fact that the classical buckling load, as well 
as that computed in the Appendix, is calculated under the assumption 
that the shell is free to expand in the radial direction during the loading 
while the edges of the test cylinders were rigidly clamped. 
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APPENDIX 


DETERMINATION OF BUCKLING LOADS FOR A CYLINDRICAL 
SHELL WITH AN AXIAL SYMMETRIC IMPERFECTION HAVING 
THE SHAPE OF A HALF SINE WAVE IN THE LENGTH DIRECTION 


The method of solution is as follows: 

1. Determine the stresses and deflections in the shell that 
occur before buckling during the loading. 

2. Consider the stresses and deflections occurring during 
the buckling as small perturbations about the solution 
of step 1 and linearize the equations. 

3. Solve the eigenvalue problem obtained in 2 and deter- 
mine the smallest eigenvalue. 

The equations that will be used for the solution are the shallow 
shell equations of Marguerre (Ref. 7). Using the coordinate system 
shown in figure 22, the equations are: 

V 4 F = Et |w 2 - w w + 2w w -w w -ww "1 (la) 

(^ xy xx yy oxy xy oxx yy xx oyyj ' 

D V7 4 w = F (w + w ) - 2F (w + w ) + F (w + w ) (lb) 

V yy xx oxx xy' xy oxy' xx yy oyy 


where w is the initial deviation from the flat plate and F is the stress 
function defined as: 
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Equation la is a compatibility relation and equation lb is the equilibrium 
equation in the direction perpendicular to the xy plane. 

In the case of interest here as shown in Figure 22: 


w 


o 




sin 


TTX 

~r 


(3) 


Let F , and w be the solutions of the axial symmetric problem before 
the buckling occurs. Let F and w be the perturbation stress function 
and radial deflection occurring during the buckling. Then 

F=F+F w = w^ + w. (4) 

* 

Substituting Equation 3 and 4 into Equation 1 and remembering that the 
state is axially symmetric, one obtains: 
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Since F and w are the solutions for the pre-buckled state, the 
terms that are underlined are equal to zero. Also, since we are only 
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interested in the bifurcation points of the solutions, the perturbations 
may be considered arbitrarily small and all nonlinear terms in F and 
w can be dropped. This leaves us with the following set of equations: 
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t rx v 

sin X” 



dV 
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F 

yy 
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F* w + F ^ + F* w , 
yy xx xx R xx yy’ 


( 6 ) 


which is seen to have a solution of the form: 


w * ^ (x) sin , F * t| (x) sin 2^ . (7) 

This solution satisfies the_necessary periodic conditions in the circum- 
ferential direction when m is a whole number. 


Axially Symmetric Solution 


The axially symmetric solution must satisfy the following set of 
equations: 



Et 
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XX 
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dV 
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_ ■ , * . ,tt ,<2 . irx. 

- F -> w xx + a o ( r> sin x> - 5 


1 f* 

R r xx 


(8b) 


It will be assumed that the shell is free to expand radially at the 
ends during the loading and that the edges of the shell are supported in 
a pinned manner as the shell expands radially. This expansion is due to 

- p R<r 

the Poisson ratio and is equal to — g — , where <r is the applied 
stress. 


11 



Therefore, the boundary conditions on 

* * R % * 

w - — ~ — e 0, w « 0 at x = 

E xx 

# 

The boundary conditions on F are given by: 


w are as follows: 
0,L (9) 
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yy 
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cr t 
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at x = 0,L 


( 10 ) 


The other conditions necessary to completely specify the problem are 

* * * 
that u and w must be independent of y and that v must be 

identically zero. 

Using these conditions the solution of the axially symmetric 
problem is given by: 
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The middle term is equal to two when <r Q is the classical buckling 

stress and the first and last terms are reciprocals of each other. The 
value of the last term, for the shells that were tested in the experimental 
part of this work, is 1832. Since the <t q of interest is always of the 

order of the classical buckling stress the denominator can be approxi- 

2 

mately rewritten as Et/R . Therefore, the quantities of interest, which 
occur in the perturbation equations, can be written as: 
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It is noticed that w occurs in the coefficient of the perturbation 

xx 2 

equations always in combination with a^ir/L) sin -j- . Therefore, this 
coefficient can be written as: 


/TT x 2 . TTX (\ °"o /RTT.2) / a . 

a 0 <i> sin ir 5 1 - it * j (14) 

For the shells of concern here, R is of the order of L and o' is 
much less than E. Therefore, this second term in the brackets can be 
neglected in comparison with one. This amounts to neglecting the 
additional curvature of the shell in the x-direction caused by the loading 
in comparison with the initial curvature of the shell in this direction. 


Solution of the Perturbation Equations 


After the initial load is applied, the ends of the shell are assumed 
to be pinned with respect to the radial direction. This implies 
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The boundary conditions that will be applied on F are as follows: 
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These conditions on F imply that at the ends of the shell the stress in 
the axial direction and v* are equal to zero. 


These boundary conditions are commonly used in shell stability 
analyses because of their simplicity. They can be looked upon as 
allowing the ends of the shell to warp freely in the axial direction, 
constraining the circumferential movement, and giving a pinned support 
with respect to the radial direction. 


Now let us introduce the nondimensional quantities as follows: 
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Using equations 13 and 17 the perturbation equations given by equation 
6 become: 
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etc. 


The solution given here will be an approximate one which satisfies 
the compatibility equation exactly and the equilibrium equation approxi- 
mately. The procedure is as follows. 


1. Assume a suitable form of £ which contains a 
number of arbitrary constants and satisfies 
exactly the necessary boundary conditions. 

2. Substitute £ into the compatibility equation and 
solve for r| making sure that r\ satisfies all the 
necessary boundary conditions. 

3. Substitute and r] obtained in the first two steps 
in the equilibrium equation and determine the un- 
known constants in such a manner so as to 
minimize the error by the Galerkin method. 

It is assumed that can be represented by the following series 
which satisfies the necessary boundary conditions: 
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Next ^ and r| are substituted in the equilibrium equation which is 
written in the following manner: 

2 rn + Xa^- a-r/+ \a jin^ sin£ + a^m^ri sin£ = E. ( 22 ) 

E represents the error resulting from the fact that this equation is not 
identically satisfied. The unknown constants are determined by a 
Galerkin procedure in the following manner: 


J E sin i £ d§ = 0, 
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( 23 ) 


The resulting equations for the A^'s then become: 


n 


/ (A. . - XB..)A. s 0 

J 

where 

O a l m ' 

B.. * Q.. - (b. . - c..) 5— 

ij ij ij 7 


i f j both even or odd, 


( 24 ) 


am 


* 


a . .2 a,m 

" + + 


mj 4 a j 




2 4 (■ _ 

a i m }C _J C 1 , r 1 ? 

4 ai 2 l bi ’ j - 2 a mj-l “i.Z-j a mj-l ° i '^ +Z a mj+l J 


, 2 r 2 

2a,m \ m e . e . J a 

&m J _Z?1L ( + 


a . m | 

TT i ( 


mi 


/%! 


mj- 
2 3 . 

Trai C 


g .h . - 2m d . g 
& mj mi mj B mi 


.} 


1 j. 1 

c. . = TT-rr-T + ' 


_ , 2 .2.2 , 1 , 1 

a mj“ ' m J ’ b ij i+j-1 i-j+1 * v 'ij “i+j+l T i-j-1 

i_.e .. f - 1 

9 m i 21 a 9 r 


d . = — - 

mi a . . a . , , ' mj a . , a . . , ' mj a . . . a . . 

J mj-1 mj+1 J mj-1 mj+1 J mj+1 mj-1 


’mj 


V ♦ n 

mj-1 


S' 


1 

1 

, h 

mj-1 

/ a mj+l 

* mj 

jL 4. Jiiii! . 

k 

a . 
mj 

a mj-l 

mj+1 


- e 


m J 


2 . 
- m d 


.2 


mj 


mj 


mr 

a mj+ 1 


17 



The minimum eigenvalue will be dependent on the number of 
circumferential waves m as well as the parameters of the problem 
a^ and a which are the amplitude of the initial imperfection and the 
properties of the shell. The solutions of the eigenvalue problem were 
found by using an IBM 7090 digital computer. 

In all the calculations carried out the parameter containing the 
shell dimensions was fixed at the value corresponding to the experi- 
mental test cylinders. This value of a was 18 32. Both the even and 
odd problems were computed. The results of this calculation for five 
values of a^ < 0 are shown in Figure 23. The classical value for 
the perfect shell corresponds to 1. The numbers on the curves 
correspond to the predominant term in the eigenvector. 

The absolute minimum eigenvalue over the whole range of the 
circumferential wave number is easily found from these curves as 
well as the corresponding circumferential wave number. The mode 
shape of this minimum is seen to be the first mode. This means that 

the buckling shape for a^ ^ 0 has the form sin in the axial 

direction. The higher modes contribute only about a 2 per cent 
correction to this pattern. The number of circumferential waves 
varies from 6 to 8 as a^ varies from -100 to 0. 

In the case where the shell bows outward, meaning aj> 0, 
convergence of the eigenvalue problem was not obtained. However, the 
classical buckling load corresponding to the value of A * 2 was 
established as a lower bound for all values of a^ 0. 

The variation of the theoretical buckling load with a Q /t is 
shown in Figure 21. 

Since for a ^ 0 the reduction in buckling load given by the theory 
is much less than tfiat found in the tests the following comparison was 
made. The value of the buckling stress was found that corresponded to 
approximately the same number of waves in the circumferential direction 
as that found in the tests. The comparison for a Q < 0 is shown in 
Figure 24. The numbers on the curves correspond to the predominant 
term in the eigenfunction expansion, and those at the ends of the curves 
to the number of circumferential waves. It is seen that the decrease in 
buckling stress with increasing negative imperfection amplitude is less 
as the circumferential wave number increases. The buckling mode shape 
is also more in line with the experimental results. The eigenfunction has 
more axial waves as the number of circumferential waves increases and 
the amplitude of the waves at the center of the shell is greater than the 
amplitude of the waves at the edges. Again, this is nearer the experi- 
mental results. 
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Fig. 1 - Mandrel and Finished Wax Form 
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STRAIN CURVE FOR PLATED COPPER 




Fig. 3 - Testing Machine 



Fig. 4 - Load Measuring Cylinder 
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Fig, 7 - Initial Imperfection Measuring Equipment 




28 








6 


8 


FIG. 


Distance Along Shell Axis (inches) 
- SHELL 7 INITIAL IMPERFECTION 




Initial Imperfection (inches) 


280 °, 320 ° 


240 ° 


\- 0.020 



120 °, 160 ° , 200 ° 




FIG. 10- SHELL X INITIAL IMPERFECTION 
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FI6. 14- SHELL AD INITIAL IMPERFECTION 
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FIG. 15- SHELL S INITIAL IMPERFECTION 
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FIG. 16 - SHELL R INITIAL IMPERFECTION 
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FIG. 19 -LOAD DISTRIBUTION NEAR BUCKLING 

SHELLS 2,3, AB, AD, S 


40 




Shell R 92% a 


Initial 

Buckle 


Shell Q 90% a 


Initial 

Buckle 


» 1.40 


Shell P 92% or 


Shell 4 95 %a 


2.20 

2.00 

1.80 

1.80 

1.60 

1.40 

1.20 

1.00 

0.80 

0.60 



Shell T 96% <j c 

Initial / 

Initial 

Buckle / 

\Buckle 


Circumferential Position in Degrees 

FIG. 20“ LOAD DISTRIBUTION NEAR BUCKLING 
SHELLS T, 4, P, Q, R 







Theory 


C I; C(E 

.2 2 . 2 ^ 

4 — O *- *• 

O > O U 
ml-® 

Z 3 **- IO 
* ° SO 
at*® 0 

6 oEo-iS 

® c — 

- O .O 
WO O 

o □ - £ 


□ 

Q> 

3 


§ 


o 

E 


o 

z 

o 

and 

V 

3 

0) 

+- 

°1§° 

00/0 

E 

a> 

c 

0 
y 

1 


pfooo 

CD 


O <» 

— w 

H UJ 

< z 

z * 

K O 

o - 


0 > N 

°s r 

w «*- 

® E « 

•O 3 ® 

E o > 

3 - 2 
z o 5 



43 


IN THE LENGTH DIRECTION 
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